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^>The  program  HEAVE  is  a  one-dimensional  computer  program  for  approximate 
analysis  of  vertical  movements  of  swelling  foundation  soils  beneath  permanent 
structures  caused  by  changes  in  vertical  loads  and/or  the  moisture  profile. 

The  program  is  applicable  to  slab,  long  continuous,  and  circular  shaft  founda¬ 
tions.  Results  of  both  one-dimensional  consolldometer  swell  tests  and  soil 
suction  tests  may  be  used  in  HEAVE., 
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20.  ABSTRACT  (Continued). 

~>>rhe  program  considers  effects  of  soil  overburden  pressures  and  struc¬ 
tural  loads,  heterogeneous  soils,  and  saturated  or  hydrostatic  equilibrium 
moisture  profiles  on  the  computed  heave.,.  An  arbitrary  final  soil  suction 
profile  may  also  be  input  if  soil  suction  results  are  available.  Differential 
heave  may  be  estimated  by  comparing  heaves  computed  for  different  soil  strata 
and  vertical  loads. 

' —  REAVE  can  estimate  the  movement  of  circular  drilled  shaft  foundations, 
the  maximum  tension  force  in  circular  shafts  from  heave  of  surrounding  swelling 
soil,  and  the  restraining  force  provided  by  the  underream  (a  bell  or  enlarge¬ 
ment  of  the  shaft  base).  Upward  movement  is  assumed  negligible  if  restraint 
exceeds  the  uplift  thrfuTt  from  surrounding  swelling  soil.  Estimates  of  shaft 
movement  can  be  made  for\  moisture  migrating  down  the  shaft  from  the  ground 
surface,  moisture  migrating  from  an  Intermediate  zone  such  as  from  a  relatively 
thin,  pervious  sandy  stratum,  and  moisture  migrating  upward  from  below  the 
base  of  the  shaft  such  asi  from  capillary  rise  of  a  rising  water  table. 
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Preface 


This  computer  program  and  user's  guide  were  prepared  under  RDT&E 
Work  Unit  AT40/E0/004,  "Foundations  on  Swelling  Soils,"  sponsored  by  the 
Office,  Chief  of  Engineers,  U.  S.  Army.  The  investigation  on  which  this 
report  is  based  was  conducted  during  FY  81. 

This  user's  guide  was  prepared  by  Dr.  I.  D.  Johnson,  Research 
Group  (RG) ,  Soil  Mechanics  Division  (SMD),  Geotechnical  Laboratory  (GL), 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  under  the  general 
supervision  of  Mr.  C.  L.  McAnear,  Chief,  SMD,  and  Dr.  W.  F.  Marcuson  III, 
Chief,  GL.  Dr.  P.  F.  Hadala,  Assistant  Chief,  GL,  Dr.  E.  B.  Perry,  RG, 
SMD,  and  Mr.  G.  B.  Mitchell,  Chief,  EG,  SMD,  reviewed  the  report. 

The  Commander  and  Director  of  WES  during  preparation  of  this  report 
was  COL  Tilford  C.  Creel,  CE.  The  Technical  Director  was  Mr.  Fred  R. 
Brown. 
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U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


_ Multiply _ 

feet 

pounds  (force)  per 
square  foot 

square  feet 

tons  (force) 

tons  (force)  per 
square  foot 

tons  (mass)  per 
cubic  foot 


_ gy 

0.3048 

47.88026 

0.09290304 

8896.444 

95.76052 

320.3692 


To  Obtain 

metres 
pascals 

square  metres 
newtons  • 
kilopascals 

kilograms  per 
cubic  metre 


USER'S  GUIDE  FOR  COMPUTER  PROGRAM  HEAVE 


Purpose  and  Scope 

1.  The  purpose  of  this  user's  guide  is  to  document  the  computer 
program  HEAVE,  which  was  developed  to  assist  the  design  engineer  in  de¬ 
termining  optimum  foundations  in  expansive  soils.  This  program  may  be 
used  to  estimate  vertical  foundation  movements  for  the  methodology  de¬ 
scribed  in  TM  5-818-7  (Headquaters ,  Department  of  the  Army  1982).  The 
program  HEAVE  approximates  potential  vertical  movements  beneath  the  cen¬ 
ter  or  corner  of  rectangular  slab  foundations,  beneath  the  central  axis 
or  edge  of  long  (or  strip)  continuous  footings,  and  beneath  the  center 
of  circular  deep  shaft  foundations.  Vertical  movement  may  also  be  com¬ 
puted  for  soils  beneath  these  foundations  constructed  in  excavations  and 
for  soils  adjacent  to  the  foundation  not  confined  by  structural  loads. 
Foundation  movement  estimated  by  HEAVE  has  been  checked  by  comparison 
with  field  test  sections  (Johnson  1981)  and  practical  application  to  ex¬ 
isting  structures. 

2.  Structural  loads  are  assumed  to  be  transferred  to  the  support¬ 
ing  soils,  and  the  soil-bearing  capacity  is  assumed  not  to  be  exceeded 
beneath  the  foundations.  Computation  of  potential  heave  beneath  the 
circular  deep  shaft  foundation  also  considers  the  restraining  force  of 
underreams  and  the  uplift  force  of  swelling  adjacent  soil.  Laboratory 
test  data  from  consolidation  swell  or  soil  suction  tests  are  used  in  the 
program  HEAVE. 


Methodology  of  Computations 


Foundation  movement 

3.  The  potential  total  vertical  heave  at  the  bottom  of  the  founda¬ 
tion  as  shown  in  Figure  1  is  determined  by 


>  FOUNDATION  SLAB 


UW,  WET  X 
PROFILE  — ^sx 

uw,  DRY  f 
PROFILE S 

METHOD  2: 

Uw  =  yUX  ■  X.) 


METHOD  1:  © 
uw  =  0  \  yi 


PORE  WATER  PRESSURE 


:  °x  \  \ 

Uw*  DRY  \ 
PROFILE 

NODAL  POINT  16 


NEL  =  21 


'  NODAL 
POINT  22 


PORE  WATER  PRESSURE 


a.  Shallow  Groundwater  Level 


b.  Deep  Groundwater  Level 


Figure  1.  Schematic  for  computation  of  vertical  movement  by 

program  HEAVE 
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AH  =  N  •  DX 


DELTA(i)  =  N  •  DX 


i=NBX 


i — NEL  / ‘ \  / . \ 

v  ef(l)  "  e0(l) 
Z  1  +  e  (i) 

i=NBX  ° 


where 


AH  =  potential  vertical  heave  at  the  bottom  of  the  founda¬ 
tion,  ft* 

N  =  fraction  of  volumetric  swell  that  occurs  as  heave  in 
the  vertical  direction 

DX  =  increment  of  depth,  ft 

NEL  =  total  number  of  elements 

NBX  =  number  of  nodal  point  at  bottom  of  the  foundation 

DELTA(i)  -  potential  volumetric  swell  of  soil  element  i  ,  fraction 

e^(i)  =  final  void  ratio  of  element  i 

e  (i)  =  initial  void  ratio  of  element  i 
o 

The  AH  is  the  potential  vertical  heave  beneath  a  flexible,  unrestrained 
foundation.  The  bottom  nodal  point  NNP  =  NEL  +  1  is  often  set  at  the 
depth  of  the  active  zone.  The  program  HEAVE  assumes  that  N  =  1  .  The 
fraction  N  should  be  1  for  one-dimensional  consolidometer  swell  test 
results.  The  soil  suction  test  results  tend  to  provide  an  upper  bound 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ments  to  metric  (SI)  units  is  presented  on  page  3. 
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estimate  of  the  maximum  in  situ  heave  for  intact  soil  (N  =  1)  in  part 
because  the  soil  suction  tests  are  performed  without  the  horizontal  re¬ 
straint  on  soil  swell  that  exists  in  the  field  and  during  consolidometer 
tests.  Dividing  the  heave  by  3  (N  =  1/3)  will  tend  toward  a  lower 
bound  estimate  for  fissured  soil  (Richards  1967).  Differential  heave 
may  be  estimated  by  comparing  the  computed  heaves  for  different  soil 
strata,  vertical  loads,  and  equilibrium  moisture  profiles. 

4.  The  program  HEAVE  computes  the  difference  in  void  ratio  Ae 

(i.e.,  e,  -  e  )  from  consolidometer  swell  test  results  of  each  soil 
f  o 

layer  by 

a 

Ae  =  c  log  -4  (2) 

s  v 

where 

c  =  swell  index 
s 

Os  =  swell  pressure,  tsf 
cT  =  vertical  effective  pressure,  tsf 
The  compression  index  c£  is  substituted  for  cg  if  exceeds  ag  . 

5.  The  program  HEAVE  computes  Ae  from  soil  suction  test  results 
of  each  soil  layer  by 


(3) 


where 


aGs/100B,  suction  index 
compressibility  factor 
specific  gravity 
slope  soil  suction  parameter 

initial  matrix  suction  without  surcharge  pressure,  tsf 
final  matrix  suction  without  surcharge  pressure,  tsf 


The  compressibility  factor  a  is  the  ratio  of  the  change  in  volume  for 
a  corresponding  change  in  water  content  and  can  be  found  by  the  pro¬ 
cedure  in  TM  5-818-7.  The  B  parameter  is  the  slope  of  the  matrix  soil 


suction-water  content  relationship  expressed  as 


log  r°  =  A  -  Bw 
®  m 


(4 


where 

A  =  ordinate  intercept  soil  suction  parameter,  tsf 
w  =  water  content,  percent  dry  weight 

6.  The  matrix  suction  t°  is  assumed  to  be  related  to  the  pore 

ID 

water  pressure  by 


t°  =  -u  +  Of  o 
m  w  m 


(5 


where 

uw  =  pore  water  pressure,  tsf 

=  total  mean  normal  confining  pressure,  tsf 
The  total  mean  normal  confining  pressure  cr^  is  given  by 

(1  +  21^) 

a  = - — - o 

m  3  v 

where 

K.J.  =  ratio  of  total  horizontal  to  vertical  stress  in  situ 
ct^  =  total  vertical  pressure,  tsf 

Equilibrium  pore  water  pressure 

7.  The  accuracy  of  the  estimates  of  the  potential  vertical  heave 
in  simulating  the  maximum  in  situ  heave  depends  heavily  on  the  ability 
to  properly  estimate  the  equilibrium  pore  water  pressure  profile.  Fig¬ 
ure  1  indicates  examples  of  extremes  that  can  occur  in  the  seasonal 
moisture  profile.  Seasonal  heave  between  extreme  wet  and  dry  moisture 
profiles  can  be  estimated  by  taking  the  difference  between  heaves  com¬ 
puted  for  both  extreme  wet  and  dry  profiles  (Figure  la),  or  the  sum  of 
the  settlement  for  the  wet  profile  and  heave  of  the  dry  profile  (Fig¬ 
ure  lb).  The  program  HEAVE  considers  the  three  following  profiles  il¬ 
lustrated  in  Figure  1. 


8.  Saturated  profile.  The  saturated  profile,  Method  1,  assumes 
that  the  in  situ  pore  water  pressure  is  zero  within  the  active  zone  X& 
of  moisture  change  and  heave 


u  =  0 
w 


9.  Hydrostatic  I.  The  hydrostatic  I  profile,  Method  2,  assumes 
that  the  pore  water  pressure  at  depth  X  becomes  more  negative  with  in¬ 
creasing  vertical  distance  above  the  groundwater  level  in  proportion  to 
the  unit  weight  of  water 


u  =  y  (X  -  X  ) 
w  *w  a 


where 


Yw  is  the  unit  weight  of  water  (0.0312  tcf). 


10.  Hydrostatic  II.  The  hydrostatic  II  profile,  Method  3,  is 
similar  to  the  hydrostatic  I  profile  except  that  a  shallow  water  table 
does  not  exist.  The  negative  pore  water  pressure  of  this  profile  also 
becomes  more  negative  with  increasing  vertical  distance  above  the  bottom 
of  the  active  zone  X  in  proportion  to  the  unit  weight  of  water 
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u=u  +Y(X-X) 
w  wa  w  a 


where  u  is  the  negative  pore  water  pressure  in  tons  per  square  foot 
wa 

at  depth  Xfl  in  feet. 

Depth  of  the  active  zone 

11.  The  depth  of  the  active  zone  Xg  is  defined  as  the  least 
soil  depth  above  which  changes  in  water  content  and  heave  occur  because 
of  climate  and  environmental  changes  after  construction  of  the  founda¬ 
tion.  The  depth  X  may  be  estimated  by  procedures  described  in 

a 

TM  5-818-7.  Predictions  of  shaft  movement  can  be  made  for  the  three 
cases  of  active  depths  shown  in  Table  1.  The  three  cases  are  differen¬ 
tiated  in  program  HEAVE  by  denoting  the  depths  Xg  and  X^  where  X^ 
is  the  depth  of  inactive  or  nonswelling  soil  overlying  the  swelling  soil. 
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Effect  of  uplift  forces 

12.  The  program  HEAVE  ignores  uplift  thrust  of  swelling  soils 
adjacent  to  slab  or  long  continuous  footings.  The  uplift  thrust  of  ad¬ 
jacent  swelling  soils  on  deep  foundations  (Figure  2)  is  determined  for 
the  assumption  that  interaction  of  stresses  between  skin  friction  and 
end  bearing  is  negligible  and  is  expressed  as 


Q  =  7ID 
xu  s 


f  dL 
s 


where 

Qu  =  maximum  uplift  thrust  on  perimeter  of  shaft,  tons 

D  =  diameter  of  shaft,  ft 
s 

L  =  thickness  of  the  swelling  layer  moving  up  relative  to  the 
n  shaft,  ft 

f  =  skin  resistance,  tsf 
s  ’ 

dL  =  differential  increment  of  shaft  length  L  ,  ft 

The  point  n  in  Figure  2  is  the  neutral  point.  The  value  of  Lq  should 
be  approximately  equal  to  the  depth  X  .  The  maximum  tension  force  T 

3 

in  the  shaft  is  estimated  from 


T  =  Qw  -  Qu 

where  is  the  loading  force  from  the  structure  and  includes  the 

weight  of  the  shaft. 

13.  The  skin  friction  f  is  evaluated  by 


f  =  c  +  Ka'  tan  d> 
s  a  v  Ta 


where 


c  =  adhesion,  tsf 
a  ' 

K  =  ratio  of  horizontal  to  vertical  effective  stress 

O'  =  vertical  effective  stress,  tsf 
v 

0  =  angle  of  friction  between  the  soil  and  shaft,  deg 


T  =  negative  tension  force,  tons 

D  =  shaft  diameter,  ft 
s  * 

15.  The  bearing  capacity  of  the  soil  above  the  underream  is  cal¬ 
culated  in  program  HEAVE  by  the  equation 


q,  =  cN  +  o'N 
c  v  q 


where 

q^  =  bearing  capacity  above  the  under ream,  tsf 
c  =  cohesion,  tsf 

=  effective  vertical  overburden  pressure,  tsf 
Nc,Nq  =  bearing  capacity  factors 

If  undrained  strength  data  are  used,  then  N  =  7  and  N  =  0  .  The 

c  ■  q 

cohesion  c  is  the  undrained  shear  strength  cu  .  If  drained  strength 
data  are  used,  then  the  bearing  capacity  factors  are  calculated  by  pro¬ 
gram  HEAVE  by  the  following  equations  (Vesic  1977): 

N  =  cot  #*(N  -  1)  (16 

c  q 

Nq  =  (1  +  tan  *,)etan  ♦'  tan2  ^45  +  (17 


where  <f>'  is  the  effective  angle  of  shearing  resistance.  The  cohesion 
c  is  the  effective  cohesion  in  Equation  15.  Equation  17  is  indicative 
of  failure  by  the  Terzaghi  local  shear  mechanism  and  tends  to  represent 
a  lower  limit  of  the  bearing  capacity  factor  Nq  . 

16.  The  soil  suction  model  described  in  TM  5-818-7  may  also  com¬ 
pute  the  effect  of  uplift  as  well  as  the  consolidometer  swell  model  if 
the  swell  pressure  is  input  for  each  soil  layer.  If  the  swell  pressure 
is  not  input,  the  program  assumes  that  the  swell  pressure  is  given  by 
the  initial  soil  matrix  suction,  an  assumption  that  can  greatly  exceed 
the  swell  pressure.  Computed  shaft  uplift  movements  will  be  an  upper 
limit  because  the  shaft  stiffness  is  assumed  equal  to  the  soil 
stiffness. 


X 

17.  The  program  HEAVE  computes  the  bearing  capacity  of  the  founda¬ 
tion  soil  assuming  that  Equations  15-17  are  applicable.  If  undrained 
strength  data  are  used,  N£  =  9  .  The  message 

BEARING  CAPACITY  OF  TSF  EXCEEDED  FOR  ELEMENT 

is  printed  if  loading  pressures  exceed  the  computed  bearing  capacity  of 
the  given  element. 


Computer  Program 

18.  The  program  consists  of  a  main  routine  and  five  subroutines. 

The  main  routine  feeds  in  the  input  data,  calculates  effective  over¬ 
burden  pressure,  determines  the  restraining  force  for  a  deep  shaft 

foundation,  and  computes  the  foundation  force  beneath  the  shaft  and  ten¬ 
sion  force  from  uplift  of  adjacent  swelling  soils.  The  subroutine  MECH 
applies  the  mechanical  model  for  prediction  of  potential  heave  using  the 
results  of  consolidometer  swell  tests.  The  subroutines  SUCT  and  HSUCT 
apply  the  soil  suction  model.  The  subroutine  PSAD  sets  up  the  proper 
depths  in  the  soil  profile  for  calculation  of  heave.  The  subroutine  SLAB 
sets  up  the  bearing  pressure  for  slab  and  long  continuous  footings.  The 
program  is  set  to  consider  up  to  10  different  soils  and  a  maximum  of 

80  soil  elements.  The  capacity  of  the  program  may  be  increased  by  ad¬ 
justing  the  PARAMETER  statements. 

Input  data 

19.  The  program  was  prepared  for  time-sharing  on  the  Honeywell 
series  G600  computer.  The  input  data  are  entered  in  free  field  format 
as  illustrated  in  Table  2.  Descriptions  of  the  input  data  are  given  in 
Table  3. 

Output  data 

20.  The  output  data  are  illustrated  in  Table  4.  Descriptions  of 
the  output  data  are  given  in  Table  5. 


Example  Applications 

21.  Three  example  problems  are  provided.  All  use  the  soil  proper¬ 
ties  given  in  Table  6. 

Slab  in  excavation 

22.  Table  7  illustrates  the  input  data  for  a  slab  in  an  excava¬ 
tion  12  ft  deep  (Figure  3).  The  groundwater  DGWT  is  22  ft  below  ground 
surface  or  10  ft  below  the  slab.  The  mechanical  model  was  used  for  a 
100-  by  100-ft  slab.  The  output  data  are  provided  in  Table  7  for  a 
loading  pressure  of  0.072  tsf  assuming  a  saturated  final  pore  pressure 
profile  at  the  center.  Heave  beneath  the  slab  is  0.25  ft.  Heave  ad¬ 
jacent  to  the  slab  in  soil  above  12  ft  of  depth  is  0.18  ft. 

Long  continuous  footing 

23.  Table  8  illustrates  the  input  data  for  an  infinitely  long 
continuous  footing  100  ft  wide  on  the  ground  surface  (Figure  4)  using 
the  mechanical  model.  The  groundwater  level  is  8  ft  below  the  footing. 
The  output  data  are  provided  in  Table  8  for  a  loading  pressure  of 
0.072  tsf  assuming  a  saturated  final  pore  pressure  profile.  The  heave 
beneath  this  footing  is  0.14  ft. 

Circular  shaft 

24.  Table  9  illustrates  the  input  data  for  a  circular  shaft  30  ft 

deep  and  2  ft  in  diameter  with  a  3-ft-diam  underream  (Figure  5)  using 
the  soil  suction  model  with  drained  strength  data.  The  active  zone  ex¬ 
tends  from  the  ground  surface  to  10  ft  below  the  base  of  the  shaft.  A 
10-ton  loading  force  is  placed  on  the  shaft.  The  output  data  are 

provided  in  Table  9  assuming  a  saturated  profile  for  the  full  40  ft  of 
active  depth,  X  The  shaft  heave  of  1.16  ft  is  the  maximum  that  a 

d 

cracked  shaft  could  rise  excluding  subsoil  movement,  while  one-third  of 
this  heave  or  0.4  ft  is  more  likely  in  a  fissured  soil.  The  actual 
shaft  heave  will  be  less  if  no  fracture  results  in  the  shaft.  The  maxi¬ 
mum  soil  heave  adjacent  to  the  shaft  is  1.3  ft.  A  gap  could  open  be¬ 
neath  the  bottom  of  the  shaft.  The  maximum  tension  force  of  -166.5  tons 
will  require  1.3  percent  steel  in  the  2-ft-diam  shaft. 
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Table  2 

Input  Data  Format 


[;r 


ft 


Ste£ 

1 


2 


3 


4A 


4B 


Data 


The  program  will  print: 

TITLE? 

=  .  Input  description  of  the  problem 

The  program  will  print  after  carriage  return: 

NPROB,NSUCT ,NBPRES ,NNP,NBX,NMAT,DX 

=  .  Input  the  above  variables  (see  Table  3  for  definitions) 

The  program  will  print  after  carriage  return: 

M,G,WC,EO,C,PHI,K 

=  .  Input  the  above  variables 

If  NSUCT  =  0  ,  the  program  will  print  after  cazriage  return: 
M,SP,CS,CC 

=  .  Input  the  above  variables, 

If  NSUCT  =  1  ,  the  program  will  print  after  carriage  return: 
M,A,B,ALPHA,KT,PI ,SP 

=  .  Input  the  above  variables.  If  ALPHA  left 

blank,  then  a  =  0  for  PI  <  5 

a  =  0.0275PI  -  0.125  for  5  <  PI  <  40 
or  =  1  for  PI  >  40 

If  SP  left  blank,  the  swell  pressure  will  be  assumed 
the  initial  matrix  soil  suction 


5 


6 


7A 

7B 


The  program  will  repeat  steps  3  and  4  until  all  soils  from  M  -  1 
to  M  =  NMAT  are  input 

The  program  will  print  after  carriage  return: 

ELEMENT, NO.  OF  SOIL 
=  .  Input  1,1 

=  .  Input  element, 2  for  elements  in  increasing  order 

for  each  increase  in  soil  type  M 
=  .  Input  NEL,NMAT  as  the  last  and  deepest  element 

for  soil  type  M  =  NMAT 

The  program  will  print  after  carriage  return  up  to  NPROB  : 
XA,XF,DGVT, I0PTI0N.N0PT 
=  .  Input  the  above  variables 

If  NBPRES  #  1,  the  program  will  print  after  carriage  return: 

Q , BLEN , BWID , LOCATION (CENTER=0 , EDGE/CORNER= 1 ) 

=  .  Input  the  above  variables 

If  NBPRES  =1,  the  program  will  print  after  carriage  return: 
PLOAD,AF,DP,DB 

=  .  Input  the  above  variables 


cs 

4A 

cc 

«A 

N 

4S 

A 

AS 

• 

AB 

ALPHA 

AB 

KT 

AB 

PI 

Al 

SP 

AB 

ELEMENT 

S 

HO.  of  Soil 

3 

ML 

3 

NNAT 

3 

XA 

* 

KP 

6 

Dcvr 

A 

IOPTION 

6 

HOPT 

A 

Q 

7A 

RUN 

7A 

■HID 

7k 

LOCATION 

7k 

PLOAD 

79 

AP 

79 

DP 

79 

DN 

79 

Nutter  of  cM«s  with  tte  mm  soil  profile.  Loads  and  foundation  dimensions  can  be  varied 

Option  for  aodel:  *0  for  consolidation  swell  (HECK)  Model;  *1  for  soil  suction  Model 

Option  for  foundation:  *1  for  circular  or  shaft;  *2  for  rectangular  slab;  =3  for  long  continuous 

Total  nunter  of  nodal  points,  NIL  ♦  1 

tenter  of  nodal  point  at  tte  button  of  tte  foundation 

Total  nunter  of  different  soil  layers 

Increnent  of  depth,  ft 


tenter  of  soil  layer 

Specific  gravity  G(  of  soil  layer  N 

Initial  water  content  wg  of  soil  layer  N  ,  percent 

Initial  void  ratio  eQ  of  soil  layer  H 

Soil  effective  cohesion  c*  or  undrained  shear  strength  c  ,  tsf 
Iff.ctiv.  «n,l»  of  Ut)[Ul  friction  CjiOlIC.c  “ 

Ratio  of  berisoatal  to  vertical  effective  pressure  K.K  *  1.0  if  left  blank 

Snell  Characterisation  by  the  Consolidation  Swell  (HECH)  Model 


Nunter  of  soil  layer 

Swell  pressure  of  soil  layer  H  ,  tsf 
Swell  index  cf  of  soil  layer  H 
Compression  index  cf  of  soil  layer  N 

Swell  Character list ion  by  the  Soil  Suction  (SUCT)  Model 

tenter  of  soil  layer 

Intercept  of  leg  suction  and  water  content  relationship  of  soil  layer  N  ,  tsf 
Slope  of  log  suction  and  water  content  relationship  of  soil  layer  H 
Compressibility  factor  0  of  soil  layer  N 

Ratio  of  total  horizontal  to  vertical  pressure  Kj  of  soil  layer  H 

Plasticity  index  PI  of  soil  layer  H  ,  percent 

Swell  pressure  a  of  soil  layer  H  ,  tsf.  O  «  i*  if  left  blank 


Nunter  of  soil  elenent 
Number  of  soil  layer  H 
Total  nunter  of  soil  elements 
Total  nunter  of  soil  layers 


Element  Characterisation 


milibriuM  Moisture  Profile 
tenant  beneath  tte  foundation  asm 


Depth  of  tte  active  zone  Eg  ,  ft.  Hotenent  beneath  tte  foundation  assumes  tte  equilibrium  moisture  profile  extends 
dew  to  tte  bottom  nodal  point  HP 

Depth  from  ground  surfeco  to  tte  depth  ttet  tte  setive  tone  begins  . »» 

Depth  to  tte  groundwater  table,  ft.  DOW  is  set  internally  t#  tte  depth  of  nodal  point  MP  if  I0PTI0R  >  l 
■quilihriun  moisture  profilet  -0  for  saturation  (Equation  7) ;  *1  for  hydrostatic  I  (Equation  I) ;  to  simulate  s 
Equation  (9).  sot  IOPTION  a  1  and  DOW  •  or  I  OPTION  e  2  if  RSUCT  *  1;  if  IOPTION  >  3  and  HSUCT  *  1,  tte 

final  total  soil  suction  without  surcharge  pressure  is  input  for  each  soil  layer 

Option  for  amount  of  output:  "•  for  forces  (if  NBPRKS  *  1)  and  total  heave;  *1  for  forces,  total  heave,  and  tte 
frectiee  and  excess  pore  pressure  at  each  depth  interval 

Usd  leg  and  Dimens  lone  ef  Fine  ill  ties 


Foundation  and  sapnrstructare  pressure,  tsf 

Radius  nf  circular  foundation  (NBftXS  ■  1);  length  ef  slab  (RBPRtt  *  2),  *0.0  for  long  continuous  (MBPRES  «  3),  ft 
0.0  for  circnlsr  foundation;  width  of  slob;  width  of  long  continuous  footing,  ft 

■0  fur  cantor  of  rectangular  alab  or  long  continuous  footing;  ■!  for  comer  of  olob  or  edge  of  long  continuous 
fasting*  **t  need  far  circular  foundation  tears  only  center  results  printed 


Lending  force  an  circular  shaft  including  weight  nf  shaft  ,  tens 
detection  factor  ef  shin  friction  tern  c(  Oquntleo  13) 

Diameter  ef  shaft  Dg  ,  ft 
Diameter  ef  hate  or  smderream  D.  .  ft 


Line 


Data 


NBPRES  j  1 

1  (If  NOPT  =  1) 

ELEMENT  DEPTH, FT  DELTA  HEAVE, FT  EXCESS  PORE  PRESSURE, TSF 

2  SOIL  HEAVE  NEXT  TO  FOUNDATION  EXCLUDING  SUBSOIL=  FEET 

3  SUBSOIL  MOVEMENT=  FEET 

4  TOTAL  SOIL  MOVEMENT=  FEET 

NBPRES  =  1 

1  FORCE  RESTRAINING  UPLIFT=  EXCESS=  TONS 

2  SHAFT  HEAVE  (EXCLUDING  SUBSOIL  MOVEMENT )=  FEET 

3  FORCE  AT  BOTTOM  OF  SHAFT=  TENSION=  TONS 

AREA  STEEL=  PERCENT 

4  (If  NOPT  =  1) 

ELEMENT  DEPTH, FT  DELTA  HEAVE, FT  EXCESS  PORE  PRESSURE, TSF 

5  SOIL  HEAVE  NEXT  TO  FOUNDATION  EXCLUDING  SUBSOIL=  FEET 

6  SUBSOIL  MOVEMENT=  FEET 

7  TOTAL  SOIL  MOVEMENTS  FEET 

8  TOTAL  SHAFT  MOVEMENTS  FEET 


Description 


ELEMENT 

DEPTH, FT 

DELTA  HEAVE, FT 

EXCESS  PORE  PRESSURE, TSF 

SOIL  HEAVE  NEXT  TO  FOUNDATION 
EXCLUDING  SUBSOIL 
SUBSOIL  MOVEMENT 

TOTAL  SOIL  MOVEMENT 


Number  of  element 
Depth  of  center  of  element,  ft 
Heave  of  increment  (element)  AH  ,  f 
Initial  negative  pore  pressure  ex¬ 
ceeding  equilibrium  pressure, tsf 
Movement  of  adjacent  soil  above  bot¬ 
tom  of  foundation,  ft 
Movement  of  soil  beneath  foundation, 
ft 

Total  heave  of  soil  adjacent  to 
foundation,  ft 


FORCE  RESTRAINING  UPLIFT 
EXCESS 


Force  restraining  uplift  ,  tons 
Qr  -  Qu  where  Qu  =  uplift  force, 


SHAFT  HEAVE  (EXCLUDING 
SUBSOIL  MOVEMENT) 

FORCE  AT  BOTTOM  OF  SHAFT 

TENSION 

AREA  STEEL 


TOTAL  SHAFT  MOVEMENT 


Movement  of  shaft  excluding  soil 
movement  beneath  shaft,  ft 
Force  exerted  on  soil  beneath  the 
shaft  (positive  quantity) ,  tons 
Maximum  tension  in  shaft  (negative 
quantity),  tons 

Percent  steel  from  Equation  14 
assuming  ASTM  A615  Grade  60  steel 
is  used.  The  minimum  steel  should 
be  1  percent  although  the  computed 
steel  may  be  less. 

Total  movement  of  shaft,  ft 


Note:  Positive  values  represent  upward  movement  or  compression  force 
while  negative  values  represent  downward  movement  or  tension  force. 
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system  tfort 

OLD  OR  NEW-OLD  HEAVE 

READY 

■♦RUN 

-SLAB  IN  EXCAVATION y  D6WT  -  22  FT  MECH 

NPROB y NSUCT y NBPRES y  NNP y NBX y NMAT y DX 

-3»0»2»45»25»4>  .5 

MyGy WCyEOyCyPHIyK 

;-.L  y  2 . 7  y  23  .  y.8yyy 

iv:  y  ALL  y  SP  y  CS  y  CC 

- i y  57  *  » 2 ♦ 2  y ♦ 045 y  .  27 

M r 0  y WC  y EO  f  C  >  PH I y  K 

••=2? 2.7V 25*  r  » 745 y  » y 

M / ALL  y SP  y  CS  t CC 

*2  r  60.  y  .c>6y  .045 »  .27 

M»G»WC»EO»C»PHI»K 

-3  y  2 . 75  y  30  .  y  .  825  y  y  y 

M y  ALL  y  SP  t CS  * CC 

”3  y 49 ♦ y . 7  y  » 03 y .27 

My  0  y WC  y  EO  y  C  y  PHI y  N 

::;4 « 2 »  V  6  y  29  ♦  y  » 828  yyy 

M  y  ALl.  y  SF'  y  CS  t  CC 

-•4  y  75  ♦  y  2  *  4  y  ♦  052  y  .  2 

ELEMENTyNO.  OF  SOIL 

•-  I  s  i 

-9  y  2 

■“17  y  3 

-25  y  4 

-44  y  4 


XA y  XF  y  DGWT  y I OPT I ON  y  NOPT 
™22. y ♦ 0  y  22 ♦ y  0  y 1 

Q y  BLEN » BWID y LOCATION  <  CENTER=0 y  EDGE/C0RNER=1 ) 
™ ♦ 072  y i 00 ♦ y 1 00 ♦ y  0 


(Continued) 


Table  7  (Concluded) 


Output  Data 


ELEMENT 

DEPTH* FT 

DELTA  HEAVE* FT 

EXCESS  PORE  PRESSURE *TSF 

1 

0.25 

0.05459 

2.18559 

2 

0.75 

0.04266 

2.15676 

3 

1.25 

0,03712 

2.12793 

4 

1.75 

0,03346 

2.09910 

5 

2.25 

0.03073 

2.07027 

6 

2.75 

0.02856 

2.04145 

7 

3.25 

0.02674 

2.01262 

8 

3.75 

0.02519 

1.98379 

9 

4.25 

0.01107 

0.41426 

10 

4.75 

0.00977 

0.38404 

11 

5.25 

0.00860 

0.35382 

12 

5.75 

0.00755 

0.32360 

13 

6.25 

0.00658 

0.29338 

14 

6.75 

0.00570 

0.26316 

15 

7.25 

0.00488 

0.23294 

16 

7.75 

0.00411 

0.20272 

17 

8.25 

0.00258 

0.21231 

18 

8.75 

0.00215 

0.18170 

19 

9.25 

0.00174 

0.15109 

20 

9.75 

0.00135 

0.12049 

21 

10  •  25 

0.00098 

0.08988 

22 

10.75 

0.00063 

0.05927 

23 

11.25 

0,00030 

0.02866 

24 

11.75 

0.00437 

0.32068 

25 

12.25 

0.04095 

2.31278 

26 

12.75 

0.03725 

2.28235 

27 

13. 25 

0.03441 

2,25191 

28 

13.75 

0.03210 

2.22146 

29 

14.25 

0.03016 

2.19101 

30 

14.75 

0.02847 

2.16054 

31 

15.25 

0.02699 

2.13005 

32 

15.75 

0,02567 

2.09955 

33 

16 . 25 

0,02448 

2.06902 

34 

16.75 

0.02338 

2.03847 

35 

17.25 

0.02238 

2.00790 

36 

17.75 

0.02145 

1.97729 

37 

18.25 

0.02059 

1.94666 

38 

18.75 

0,01978 

1.91599 

39 

19,25 

0.01902 

1.88528 

40 

19.75 

0.01830 

1.85454 

41 

20.25 

0.01763 

1,82375 

42 

20 . 75 

0.01698 

1.79292 

43 

21.25 

0.01637 

1.76205 

44 

21.75 

0.01578 

1.73113 

SOIL  HEAVE* FT*  ADJACENT  TO  FOUNDATION*  0,17570  SUBSOIL*  0.24608 


Table  8 


it  and  Output  Data  for  Lona  Continuous 


-LONG  CONTINUOUS  ON  GROUND  SURFACE,  DGWT  =  8  FT  MECH 

NPROB , NSUCT , NBPRES » NNP , NBX , NMAT , DX 

~2,0,3,17*1,2, ,5 

M»G»WC,EO,C*PHI,K 

-1 ,2,7,23* , .8, , , 

M , ALL , SP , CS , CC 
=1,57., 2. 2, ,045,. 27 
M»G»WC,EO,C*PHI»K 
--2 ,2*7,25,  ,  ,745,  ,  , 

M , ALL , SP , CS , CC 
~2,60, , .66, ,045, ,27 
ELEMENT, NO.  OF  SOIL 
—  1,1 
-9,2 
™16,2 

XA , XF , DGWT , IOPTION , NOPT 

-8. , .0,8. ,0,1 


3 , BLEN , BWID , LOCATION  <  CENTER=0, EDGE/CORNER® 1 ) 
•••-♦072,  .0,100.,! 


aut  Data 


.EMENT 

DEPTH, FT 

DELTA  HEAVE, FT 

EXCESS  PORE  PRESSURE 

1 

0.25 

0.03832 

2.13549 

2 

0.75 

0.03432 

2.10679 

3 

1,25 

0,03141 

2.07809 

4 

1.75 

0.02911 

2.04939 

5 

2.25 

0.02722 

2,02068 

6 

2.75 

0.02561 

1.99198 

7 

3. '25 

0.02420 

1,96328 

8 

3.75 

0.02296 

1.93458 

9 

4  ♦  25 

0.00903 

0.36518 

10 

4.75 

0.00794 

0.33508 

1 1 

5.25 

0.00694 

0.30498 

12 

5.75 

0.00603 

0.27489 

13 

6.25 

0.00519 

0.24479 

14 

6,75 

0.00441 

0.21469 

15 

7.25 

0,00367 

0.18459 

16 

7.75 

0.00299 

0.15450 

SOIL  HEAVE, FT!  ADJACENT  TO  FOUNDATION®  0,  SUBSOIL®  0,13968 

♦BYE 

♦♦resources  used  $  2.83,  used  to  date  *  589,92=  29% 

♦♦time  sharing  off  at  10,481  on  06/04/80 


Table  9 

Shaft  Foundation  30  Feet  Deei 


Input  Data 

RUN 

-SHAFT  FOUNDATION*  DGWT  -  40  FT  SUCT  DRAINED  DATA 
NPRQB *  NSUCT *  NBPRES *  NNP * NBX * NMAT *  DX 
=11*1*1*41*31*5*1. 


M  *  G  *  UC *  EO  *  C 

»PHI *K 

=1*2.7*23.* 

» 8  * . 0  *  20 . *  1 ♦ 

M* A»B» ALPHA 

* AKO*PI 

=1*6.75* .25 

*1. ,1. *39. 

M*G»WC*EO*C 

*PHI*K 

=2*2.7*25. » 

.745* .0*20. *1. 

M*A*B> ALPHA 

» AKO  *PI 

=2*6.75* .25 

*1**1. *40. 

M  *  G  *  UC  *  EO  *  C 

*PHI*K 

=3*2.75*30. 

*.825*. 0*30. *1. 

M» A*B* ALPHA 

* AKO* PI 

=3*5.04* .17 

*. 26*1. ,14. 

M*G*WC*EO*C 

*PHI*K 

=4*2.76*29. 

* .828*1**10. *2. 

M*A*B» ALPHA 

» AKO  *  PI 

=4*5.86* .18 

*1.*2.*55. 

M  *  6*  UC  *  EO*  C 

*PHI*K 

=5*2.76*29. 

*.828*1. *10. *2. 

M* A*B» ALPHA 

* AKO* PI 

=5*6.14* .19 

*1«*2.*55. 

ELEMENT* NO. 

OF  SOIL 

=1*1 

=5*2 

=9*3 

=  13*4 

=31*5 

=40,5 

XA*XF*DGWT* 

IOPTION*NOPT 

=40. » .0*40. 

*0*1 

PLOAD* AF,DP 

*DB 

=10»*1.*2«* 

3. 

(Continued) 


Table  9  (Concluded) 


FORCE  RESTRAINING  UPLIFT-  55.758  EXCESS*  -120.696  TONS 
SHAFT  HEAVE*  1.16250  FEET 

FORCE  AT  BOTTOM  OF  SHAFT*  -166.454  TENSION*  -166.454  TONS 


ELEMENT 

DEPTH » FT 

DELTA  HEAVE. FT 

EXCESS  PORE  PRESSURE 

1 

0.50 

0.15241 

9.97117 

2 

1.50 

0.12378 

9.91352 

3 

2.50 

0.11047 

9*85586 

4 

3.50 

0.10171 

9.79820 

5 

4.50 

0.06707 

2.90143 

6 

5.50 

0.06146 

2.84099 

7 

6.50 

0.05683 

2.78055 

8 

7.50 

0.05288 

2.72011 

9 

8.50 

0*01048 

0.56546 

10 

9.50 

0.01783 

0.72427 

11 

10.50 

0.01680 

0.70835 

12 

11.50 

0.01586 

0.69244 

13 

12.50 

0.06457 

3.62349 

14 

13.50 

0.04282 

3.01758 

15 

14.50 

0.04017 

2.91613 

16 

15.50 

0.03771 

2.81469 

17 

16.50 

0*03540 

2.71325 

18 

17.50 

0.03323 

2.61181 

19 

18.50 

0*03118 

2.51036 

20 

19.50 

0.02924 

2.40892 

21 

20.50 

0.02740 

2.30748 

22 

21.50 

0.02564 

2.20603 

23 

22.50 

0.02397 

2.10459 

24 

23.50 

0.02237 

2.00315 

25 

24.50 

0.02084 

1.90171 

26 

25.50 

0.01937 

1.80026 

27 

26.50 

0.01796 

1.69882 

28 

27.50 

0.01660 

1*59738 

29 

28*50 

0.01529 

1.49593 

30 

29.50 

0*01402 

1.39449 

31 

30.50 

0.01133 

1*19369 

32 

31.50 

0.01021 

1.09224 

33 

32.50 

0.00912 

0.99080 

34 

33*50 

0.00807 

0.88936 

35 

34.50 

0.00705 

0.78791 

36 

35.50 

0.00606 

0.68647 

37 

36*50 

0.00509 

0.58503 

38 

37.50 

0.00415 

0.48359 

39 

38.50 

0.00324 

0.38214 

40 

39.50 

0.00235 

0.28070 

SOIL  HEAVE* FT  1  ADJACENT  TO  FOUNDATION*  1.30535  SUBSOIL*  0*06666 


APPENDIX  A;  PROGRAM  LISTING 


lOOOC 

TOTOt 

1020C 

1030 

T0"40” 

10501 

10601 

TORO 

1085 

1086 
TOW 
1100 
1110 
1120 
1130 
1140 
IT50 
1160 
1170 
II  BO 
1190 
1200 
1210 
1215 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 


PREDICTION  OF  HEAVE  -  HEAVE 

RACED  ON  CONSTANT  VOLUME  SUELL/8UELL  OVERBURDEN/SUCYlON" 
DEVELOPED  BY  L.  D.  JOHNSON 
PARAMETER  NL<*10»NQ*81 

COMMON “  A ( NL  >  t B  <  NL ) »  G  <  NL  > » HC  <  NL  > » EO  <  NL ) » SP  <  NL  > » AK  < NL ) »  CS ( NL ) . 
CC(NL) » ALPHA <NL> r AKO(NL) »P(NQ) . IE<NQi 1 ) fNI »N2.NBX. NEL » IOPTION 
NOPT  *  DX » DXX >  DGWT » PI I  *  XA » XF » DELH » DELHI t DELH2 
“  DIMENSION  PP(NO) »PI (NL) *C(NL> »PHI TNL) 

PRINT  1 

1  FORMAT (6HTITLE? ) 

DEAD  3  -  *  - 

3  FORMAT (30H 
GAW=0. 03125 
PII-3< 14159265 
NP=1 
PRINT  5 

5  FORMAT (34HNPR0B»NSUCT»NBPRES»NNP»NBX»NMAT»DX) 

READ  »NPROB » NSUCT » NBPRES  » NNP » NBX  » NMAT r DX 
NEL=NNP- 1 
14  PRINT  10 

10  FORMAT ( 17HM  »G» WC»EO»  C  tPHI »K ) 

READ»M»G<M> »MC(M) > EO(M) t C(M> >PHI <M> >  AK  <  M  > 

PHI ( M)=PHI (M) *P1 1/180 .  '  "* 

IF(AK(M).LT«0«01)AK(M>*1.0 
IF (NSUCT . EQ> 1 )G0  TO  25 
PRINT  12 

12  FORMAT ( lOHMf SP>CS»CC) 

READ»M*SP(M> rCS(M)*CC(M) 

GO  TO  20 

25  PRINT  8 

8  FORMAT  <  20HM »  A , B  f ALPHA » KT . PI , SP  > 

READ* M>A<M) *B(M) r ALPHA (M) > AKO(M) rPKM)fSP(M) 

IF ( ALPHA (M> . LE . 0. )G0  TO  16 
GO  TO  20 

16  ALPHA(M) » . 02754PI <M>- . 125 
IF (PI  (M) .LE.5. >  ALPHA ( M ) *0 . 0 

IF(PI(M> .GE.40. >ALPHA<H)-1. 

20  IF (NMAT-M)24»27» 14 

26  PRINT  17»  M 

17  FORMAT  <  20H  ERROR  IN  MATERIAL  »I5> 

STOP 

27  L=0 
PRINT  30 

30  FORMAT  < 19HELEMENT >N0<  OF  SOIL) 

40  READtNt IE(Nf 1 ) 

50  L*L+1 

IF(N-L)60f60*70 
70  IE<Lf1)=IE<L-1»1> 

GO  TO  50 

60  IF(NEL-L)80tB0i40 
80  CONTINUE 
PRINT  81 


«1  fTOMAT</-.*3HXA.1<f  <©SWT»X0PT10M.HePf> 

1510  READ.XA.XF . DOUT . IOPTION.NOPT 

1515  IF<NSUCT.EQ.0.AND.I0PTI0N.GT.1)PRINT  85 

1517  85  FORMAT (30HIQPTX0N  TOO  LARGE)  XOPTION  -  ?> 

tstr-  --  rr rugucT.gq  <t)-.amp  .  topt  tow .  bt<i tread >  iopt  i  oh - 

1519  IF<I0PTI0N<0T<1 > DOUT “FLOAT (NEL)tDX 

1520  IF<NBPRES<EQ<1 )00  TO  89 

1530  PRINT  86 

1540  86  F0RHAT</.44HQ.BL£N.DMID.L0CATI0N<CENTER«0.EDGE/C0RNER*1)> 

1550  ‘  READ<Q<BLEN<BWID<MRECT 

1560  80  TO  91 

1570  89  PRINT  90 

1580  90  F0RMAT</ <14HPL0AD<AF< DP <DB)  ' 

1590  READ<PLOAD<AF<DP<DB 

1592  1F<DP.8T<DB)PR1NT  88 

1593  88  F0RMAT<51HDIAH  SHAFT  TOO  LARGE)  DIAM  SHAFT  -  ?<  DIAN  BASE 

1594  IF<DP.8T . DB)READ<DP<DB 

1600  91  IF< IOPTION.LT. 3)80  TO  92 

1610  PRINT  101 

1620  101  F0RHAT<31HMATERIAL  NO < FINAL  TOTAL  SUCTION) 

1630  102  READ.H<CC<H) 

1640  IF<NMAT-M) 103<92< 102 

1650  103  PRINT  105<M 

1660  105  FORMAT < 20H  ERROR  IN  MATERIAL  <I5) 

1670  GO  TO  102 

1690C  CALCULATION  OF  EFFECTIVE  OVERBURDEN  PRESSURE 

1700  92  P<1>*0.0 

1710  PP<1>*0.0 

1720  DXX-DX 

1730  DO  100  I~2<NNP 

1740  HTYP*IE< 1-1 < 1 ) 

1750  WCC*UC<MTYP ) /100 . 

1760  GAHH-G<HTYP)*GAW*<l.+WCC)/<i.4E0<MTYP) >  _ 

1770  IF<DXX<GT<  DGWT ) GAHH-G AMM-GAW 

1780  P<I)*P<I-1) +DX4GAMM 

1790  PP<I)-P<I)  _  _ 

1800  DXX*DXX+DX 

1810  100  CONTINUE 

1820  IF <NSUCT <GT <O.AND, I OP T I ON. EQ. 2) GO  TO  111 

1830  GO  TO  112 

1840  111  NATNEL*IE<NEL<1) 

1845  IF <  AKO < HATNEL ) . LT . 0 . 01 )  AKO <  MATNEL  > -AK  < MATNEL  > 

1850  F= < 1, +2. *AKO< MATNEL)) 73. 

1860  TFI *A < MATNEL )-B< MATNEL )$WC< MATNEL) 

1870  TFI=10.**TFI 

1872  ALP= ALPHA  <  MATNEL ) 

1874  IF  < DGWT . LT . FLOAT  <  NEL  >  *DX  > ALP« 1 . 0 

I860  TFI*TFI-P<NNP)*F*ALPHA<MATNEL>  _  _ 

1890  112  DXX-0.0 

1900  DO  114  1*1. NNP 

1910  A I = I - 1 

1920  BN*DGWT/DX-AI 


1930 

19403 


IF<N8UCT.EQ.O.AND.DXX.LT.DOWT.AND.IOPTXON.EQ.1>P<I>*P(IH’ 
BNBDXtBAW 

1730 - IFniB0rr.EQ70>OT  T0  T13 - 

1953  IF(X.EQ.l>MTYP-IE<lfl> 

1957  XF(I.8T.l>MTYP>XE<l-lfl> 

T9 70 - XLP“ALPH1TWTYP) - -  - - 

1975  XF(DXX.BT.DSNT>ALP-1.0 

1980  TF-0.0 

T9V0 - IFTTOFTrOHYEUTI .  AW.  DXX7CT .  BBWTVTFKIHtDXBlJAO - 

2000  I F  < I0PT ION . EQ . 2 ) TF-TF I  ♦ <  FLOAT  <  NEL  > tDX-DXX ) *0 AW 

2005  IF  <  AKO<  NT YP  > . LT . 0 . 01 >  AKO  <  NT YP  >  »AK  <  NT YP  > 

-ZOTO - F=n  7f2  f*  AKO  fNTYP  >)  /  3 . - 

2020  IFUQPTI0M.EQ.3>TF«CC<MTYP>-P<X>*F*ALP 

2040  P< I )*TF+P< I )BFBALP 


2030 - IFTPTT)  .XT  VO  .07 PRINT  lTg.TTOTI - 

2040  114  F0RHAT(31HN£6ATIVE  FINAL  EFFECTIVE  STRESS. F10. 5. 

20703  12H  IN  ELEMENT » 15) 

2080  113  DXX-DXX+DX 

*2090 — 114  CONI  I  HUE. - 

2100  IF < NBPRE8 . 6T . 1 ) CALL  SLAB ( 8 . NSUCT . BLEN . BWID. HRECT . NBPRE8 »PP(NBX> ) 

2110  IF(NBPRE8.8T .1)80  TO  210 

2120C  CALCULATION  OF  RESTRAINING  FORCE 
2130  CON»DX*PII*DP*AF 

•2TW - PTK07O - - 

2150  PR1*0.0 

2140  PS1-0.0 

2170 - AWKXA/BX - - 

2180  N1«IFIX(AN1>+1 

2190  N2-NBX-1 

2200 - iFna  .Br.wTOOTO  122  ™  '  - - 

2210  DO  120  I-N1.N2 

2220  HTYP-IEU.l) 

7Z30 - PH»PHICNTYP>  '  ~  "  “ 

2240  TA-8IN<PH)/C08(PH) 

2250  IF(N8UCT.EQ.0. OR. SP(NTYP).6T. 0.01)80  TO  115 

22SO - TACOMA  OTTYPT-B  (WrYlP)  *NC  ( NTYP") 

2270  8P(HTYP)-10.**TAUI 

2280  115  PS1»PS1*8P<NTYP)*C0N 

■2270  ~PR» ( pp  ( ITfFP  X  i  *"HT727 

2300  PR1«PR14PR*C0N 

2310  Pl-Pl  +  <  PRBTABAK ( NTYP ) +C (NTYP ) ) *CON 

2320  —  120  CONTINUE  - - 

2330  122  HAT-XE(NBX-l.l) 

2340  QBU-7.*C(NAT> 

Z33<T  -  IFTPHrOIADTCT . OTOOlVOO  10125 - 

2340  PH-PHI(MAT) 

2370  TA-8IN(PH)/C0S(PH) 

2380" -  XNQ»TPH72.)445Y*PTI/180.  - - 

2390  XNQ«(1.4TA)*EXP(TA)*(8IN(XNQ>/C08(XN0>)*«2. 

2400  XNC-(XN0-1. )/TA 

2410' - -  O10"C  (ffAT)  *XN(r+P('NBX)*XNG - 

2420  125  PRE-PL0AD4P14QBU*PII*<DB**2.-DP**2.)/4. 


2430 


PRF-PRE / ( <  FLOAT <  M2  >  4DX-XF ) *P 1 1 «DP  > 


J44DC 

2430 

2440 

jvrcr 

2480 

2490 

330TT- 

2310 

2320 

2330- 


2340 

2330 


2370 

2373 

7380 

2390 

2400 


CALCULATION  OF  EXCESS  RESTRAINING  THREE  ¥T  SOTTOH  OF  TIER 
P2-0.0 
DELAV-0.0 

- PHRT^OTO - 

PR3-0.0 

PS2-0.0 

- CTLL  pyup - 

DO  ISO  I-N1»N2 
HTYP-IE<I»1> 

- PHBPHXTHTTP7 - 

TA-8IM<PH)/C08(PH> 

IF(N8UCT.EQ.0>00  TO  143 

- C8TIITTPT«MRrPKATirrTPT>0THTYPl  /7T00  niTHTYFn 

TAUI-A<HTYP)-D<HTYP)*MC<HTYP> 

IF<8P<NTYP).GT. 0.01)80  TO  143 

- 8PT  HTTPT-TOTTtT  ATTI - 

143  PS2-P82*SP< HTTP) SCON 

PR-(PP< I )9PP< IF1 ) )/2. 


2420  PR2-PR2+PRSC0N 

2430  CT-CG(HTYP) 

7490 - rF(PRR.~BT.  8P  (  ltTTPr.imD  .H8UCT .  EO  70  >  CT-CC  THTTP  V 

2430  CT-CT/< 1 . +EO<NTYP} ) 

2440  I F  <  PRR . LT . PRF ) PRR-PRF 


7470  -  ’CA-BFOITTPT7PlRK - 

2480  DEL-CT«AL0G10<CA)*DX 

2490  if  <  DEL .  LT .  0 . 0 .  AND .  DXX .  OT .  DGWT .  AND .  N8UCT .  BT  .  0 )  DEL-DEL/ ALPHA  < HTTP-) 

7700 - irrPECTLT .  OVO  ;HND  7T0FTT0N .  LT .  2 .  AND .  NSUCT.  ST .  0  I'DEL-BEL/ALPHMHTYPT 

2710  P2-P2+  <  PR4TA4AK  <  HT  YP )  +C  <  NT  YP  )  > SCON 

2720  DELAV-DELAV+DEL 

2730  150  CONTINUE 

7740 - P8T-FS1+PS2 - 


2750  PRT-PR1+PR2 

2740  CAT-P1+P2 

2770  DP8R-P82-PR2 

2780  Q-PRE-DPSR 

7790 - IFTOPBir7BTn*2TOiPRr-T7 - 

2800  PRINT  140.PRE.Q 

2810  140  FORMAT <  / » 23HF0RCE  RESTRAINING  UPLIFT-.F10.3.VH  EXCESS-. 

70201 - FT0.3.6H  T0N81 

2830  IF(Q.GT . 0.0) DEL A V- 0.0 


2840  PRINT  142.DELAV 

7830"  T42  TORMAT T 4YHSHA F  TH  EA  U  E  <  EXCLUDING  "5088011"  HOVENEnTT-  . 
28321  F13.S.4H  FEET) 

2840C  CALCULATION  OF  FOUNDATION  PRESSURE  DENEATH  FOOTING 

7870—  -  PGTP1RT-"P5T=PRT  "  - - - 

2880  0Q-PL0AD-P8TPRT 

2890  IF  < PSTPRT . GT . CAT  >  OQ-PLOAD-CAT 

7893 - 1  f  com  it  ro.oTS  cr«o  .0  0000 

2900  T-PL0AD-DP8R 

2910  IF<DP8R.GT,P2)T-PL0AD-P2 


■i 


■ 


■2VT5 -  ~ir<T.BT.«.0rr«0VW00'0 

2920  PRINT  170*QQ*T 

2930  170  FORMAT < 2SHF0RCE  AT  BOTTOM  OF  SHAFT**F10.3*9H  TENSION*  * 

29401  F10.3*4N  TONS) 

2942  IF<T.LT.-0.0001)A8TEEL*-0.03*T/DP*»2. 

2944  IF(T.LT .-0.0001) PRINT  180* A8TEEL  _  _ 

2944'  180  F0RHAT<34X> UHAREA  8TEEL**F10»3*9H  PERCENT)' 

29S0  IF<QQ.LE. 0.01)60  TO  290 

2940  BPRES*QQ+PP ( NBX  >  * . 78544 <  DBBB2 . -DPBB2 « > 

2970  BPRES*BPRE8/ ( . 7854BDB882 . ) -PP  <  NBX ) 

2975  BPRE1*BPRES 

2980  DXX*0.0  _ 

2990  00  200  I*NBX*NNP 

3000  IF(I.EQ.NBX)00  TO  201 

3005  MTYP*IE<I-1.1) 

3004  IF<  NSUCT . EQ • 1 ) ALP* ALPHA ( HTTP ) 

3007  I F  <  NSUCT . EQ . 1 . ANO . DXX . OT . DQHT  > ALP* 1 . 0 

3008  IF  < NSUCT . EQ . 1 >  BPRE8-BPRE ItALPBt 1 . +2 . *AKO < MTYP ) )/3. 

3010  TP-1.* (DB/(2.*DXX>>*«2. 

3020  TP*TP*«1.5 

3030  P ( I ) *P( I ) +BPRES4 ( l.-l./TP) 

3040  60  TO  205 

3050  201  P(I)*P(I) +BPRES 

3040  205  DXX-DXX+DX 

3061  200  CONTINUE 

3062C  CHECK  FOR  BEARING  CAPACITY 

3063  210  DO  250  I*NBX*NEL 

3044  MTYP*IE(I»1> 

3065  QBU*9.»C(HTYP) 

3066  IF(PHI(HTYP).LT. 0.001)60  TO  270 

3067  PH*PHI<MTYP) 

3068  TA*SIN(PH)/COS(PH) 

3069  XNQ=<PH/2.)*45.*PII/180. 

3070  XNQ* ( 1 . +TA ) *EXP  <  T A ) * ( SIN ( XNQ  > /COS (XNQ) ) **2 . 

3071  XNC*(XNQ-1. )/TA 

3072  QBU*C ( MTYP ) BXNC+  <  <  PP  < I ) *PP  < I ♦ 1 > ) /2 . > *XNQ 

3073  270  IF<P<I) .GT.QBU)PRINT  280*QBU*I 

3074  280  FORMAT < 19HBEARING  CAPACITY  OF  * F10 .3* 18H  TSF  EXCEEDED  FOR* 

30751  8H  ELEMENT* 15) 

3076  250  CONTINUE 

3077C  CALCULATION  OF  MOVEMENT  FROM  MODELS 
3080  290  IF(NOPT . EQ.0)60  TO  300 

1100  PRINT  305  . 

3110  305  FORMAT (/*33HELEMENT  DEPTH.FT  DELTA  HEAVErFT* 

31201  26H  EXCESS  PORE  PRESSURE* TSF ) 

3130  300  IF ( NSUCT. EQ.O) CALL  MECH 

irfO  IF CNSUCT .'GT .0)CALL  SUCT  "  '  ' 

3141  DEL1-DELH1+DELH2 

3143  PRINT  306 * DELH 1  * DELH2 » DEL  1 

3144  306  FORMAT ( / * 48H80IL  HEAVE  NEXT  TO  FOUNDATION  EXCLUDING  SUBSOIL** 

31451  F8»5*6H  FEET  * / * 17HSUBS0IL  MOVEMENT*. 31X.F8. 5. 6H  FEET* 

31461  /.20HT0TAL  SOIL  H0VEMENT**28X*F8.5*6H  FEET) 


•  •  '  •'  •  ' 


A5 


3147  IF  < NBPRES . EO . 1 ) DEL2-DELAV+DELH2 

3148  IF ( NBPRES <EQ« 1 ) PRINT  308.DEL2 

3147  308  FORMAT ( 21 HTOTAL  SHAFT  H0UEM£NT->27X»F8.5'6H  FEET) 

3150  NP-NP41 

3160  IF(NP.6T.NPR0B)G0  TO  310 

3170  ’  GIT  TO  80  . 

3180  310  STOP 

3190  END 

"32000  -  ~  . .  '  ~ 

3210C 

3220  SUBROUTINE  MECH 

-3220 - PARAMETER  ML-10* NO-81 

3240  COMMON  A(NL) >B(NL) »G(NL> >UC<NL) *EO(NL> »SP(NL) * AK(NL) »C8<NL) . 
3250S  CC  <  NL ) »  ALPHA ( NL )  »  AKO  <  NL ) »P<NQ) »IE(NQ»1) »N1 »N2t NBXiNEL. IOPTION 
32301 — N0PTtOlvm0riD0WT»T'Tr»lA » XF  »DELH » DELHI  ."DELW2 
3280  DELHI-0.0 

3290  CALL  PSAD 

3300  -  *  IF(N1.0E.N2>00  TO  50 

3310  DO  10  I-N1>N2 

3320  NTYP-IE( I » 1 > 

3330  PR-(P( I )+P( 141 ) )/2 • 

3340  CA«SP(NTYP)/PR 

3350  E-E0<HTYP>4CC<MTYP>*AL0Q10(CA> 

3360  IF(PR«LT . SP(MTYP) )E*EO(MTYP  >4CS(HTYP)*AL0G10(CA  ) 

3370  DEL* <E-EO< HTTP) )/< 1 . 4E0(NTYP> ) 

3380  IF (NOPT • EG >0)60  TO  40 

3390  DELP-SP(MTYP)-PR 

3400  PRINT  200 f I > DXX > DEL r DELP 

3410  40  DELH1-DELH14DX4DEL 

3420  DXX-DXX4DX 

3430  10  CONTINUE 

3450  50  DELH2-0.0 

3470  IF(NBX.GT.NEL)GO  TO  175 

3480  DXX-FL0AT<NBX>*DX-DX/2. 

3490  DO  100  I-NBX'NEL 

3500  MTYP-IE(I.l) 

3510  PR-<P<I)4P(I41>)/2. 

3520  CA-SP(HTYP)/PR 

3530  E-EO ( MTYP ) +CC <  MTYP ) BAL0G10  <  CA ) 

3540  IF(PR.LE»SP(HTYP) ) E-EO < MTYP )4CS< MTYP )*AL0G10(CA) 

3550  DEL- (E-E0(MTYP)>/<1.4E0< MTYP)) 

3560  IF (NOPT •EQ*0)G0  TO  125 

3570  DELP«SP(MTYP)-PR 

3580  PRINT  200 > ItDXX»DEL»DELP 

3590  125  DELH2-DELH24DXBDEL 

3600  DXX-DXX4DX 

3610  100  CONTINUE 

3630  200  FORMAT (I5»F10«2»F15«5»5X*F15(5> 

3660  175  RETURN 

3670  END 

3680C 
3690C 


3700  SUBROUTINE  SUCT 

3710  PARAMETER  NL*10tNQ=81 

3720  COMMON  A<NL>  fB<NL>  »G<NL> »WC<NL> »EO(NL>  »SP(NL) » AK<NL) »CS(NL) . 
37301  CC(NL) » ALPHA (NL) » AKOCNL) »P(NQ > » IE(NQ f 1 > fNI »N2 . NBXfNEL t I0PTI6N. 
37401'  NOPT fDXfDXXfDGUTf PH fXAfXFfDELHf' DELHI fDELH2 
3760  DELHI =0 . 0 

3770  CALL  PSAD 

3780  IF(N1.GE.N2)G0  TO  50 

"3770  '  CALL  HSUCT  - 

3800  DELH13DELH 

3810  SO  DELH2=0  .  0 
3820  IF(NBX.GT.NEL)GO  TO  17S 

3830  DXX=FLOAT (NBXXDX-DX/2  • 

*3840 .  N1=NBX 

3850  N2SNEL 

3860  CALL  HSUCT 

3870  DELH2=DELB 

3910  175  RETURN 

3920  END 

3930C  _  - 

3940C 

3950  SUBROUTINE  HSUCT 

3960  PARAMETER  NL=10fNQ*81 

3970  COMMON  A(NL) fB(NL)  fG(NL) »UC<NL) fEO(NL) fSP(NL) . AK(NL) fCS(NL) f 
39801  CC ( NL )» ALPHA (NL) i AKO(NL) .P(NQ) f IE ( NQf 1 > fNI fN2fNBXfNELf IOPT IONf 
39901  NOPT fDXfDXXfDGUTfPII fXAfXF fDELHf DELHI fDELH2 
4010 

4020  DELH=0.0  _ 

'4030  DO  10  I~N1»N2 

4040  MTYP=IE< I  f  1 ) 

4045  IF  <  AKO( MTYP ) .LT . 0.01 )AKO( MTYP) * AM MTYP ) 

4050  F=a.+2.*AK0<MTYP))/3. 

4060  PR=  <P(I)+P(I+1) )/2  • 

4070  TAUI  *A  <  MTYP ) -B<  MTYP  ):1MC(  MTYP)  _ 

4080  TAUI=10.**TAUI 
4090  UINIT=TAUI-PR 

4100  CT “ALPHA ( MTYP )  *G <  MTYP )/<  100 . '.IB ( MTYP > ) 

4119  CT*CT/(1.+E0<MTYP) ) 

4120  RTAU=TAUI/PR 

4130  DEL-CT«AL0G10(RTAU)*DX 

4140  IF<DEL.LT .0.0. AND.DXX.GT . DDHT ) DEL =DEL/ ALPHA (MTYP) 

4150  IF < DEL. LT. 0.0. AND. IOPT ION. LT,2)DEL»DEL/ALPHA< MTYP) 

4160  IF (NOPT ,EQ.O)GO  TO  33 

4170  PRINT  30fIfDXXfDELfUINIT 

4180  30  FORMAT (I5fF10.2fF15.5f5XfF15.5) 

4190  33  DELH*DELH+DEL 

4200  DXX3DXX+DX 

4210  10  CONTINUE 

4220  RETURN 

4230  END 

4240C 

4250C 


A7 


4260 
4270 
4280 
42901 
43001 
4320 
'4330 
4340 
4350 
4360 
4380 
4390 
4400 
4410 
4420 
44  JOS 
4440  S 
4450 
4460  * 
4470 
4480& 
4490S 
-4310C 
4520 
4530 
4540 
4550 
'4552  ' 
4560 
4565 
“4367 

4568 

4569 
4370 
4571 
4580 
"4390 
4600 
4610 
4620" 
4630 
4640 
'4630  - 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 


SUBROUTINE  PSAD 
PARAMETER  NL=10»NQ=81 

COMMON  A(NL) iB(NL) »G(NL) »WC(NL> »EO<NL>  »SP<NL>  <AK(NL> »CS(NL)  » 

CC ( NL ) *  ALPHA <NL ) >  AKO (NL)rP(NQ)rIE(NQrl)>Nl>N2iNBX>NELf I OPT  ION  * 
NOPT  »  DX  »  DXX  »  D6HT » PI I r  XA  »  XF  »  DELH , DELHI . DELH2 
AN1=XF/DX 
AN2=XA/DX 
N1=IFIX( AN1 ) +1 
N2-AN2 

BXX=XF+DX/2. 

N3=NBX-1 

IF (N2.GT . N3)N2=N3 

CONTINUE  "  ' 

RETURN 

END 


SUBROUTINE  SLAB < Q , NSUCT » BLEN.BMID»MRECT»NBPRES»WT) 

PARAMETER  HL'=10  .HQ =81  '  -  - 

COMMON  A(NL) rB(NL) . G<NL) . WC(NL) . EQ<NL) »SP<NL) * AK(NL) »CS(NL) » 
CC(NL) »ALPHA(NL) r AKO(NL) »P(NQ) »IE(NQ» 1 ) »NltN2fNBX*NEL. IQPTIQNt 
NOPT  .  DX  r  DXX .  BGUT  <PI  I  .  XA.  XF  .  DELH  >  DELHI  r  DELH2 

CAT-CGLATION  OF  SURCHARGE  PRESSURE  FROM-  STRUCTURE - 

NNP=NEL+1 

ANBX  =  FLOAT <  NBX ) *BX 
DXX  =  0.0 
BPRE1=Q-UT 

BPRES=BPREl  ~  .  : 

DO  10  I=NBX » NNP 
IF<DXX<LT<0<0D GO  TO  30 

MTYP=IE<  1-1 » 1 )  - - 

IF ( NSUCT .  EQ .  D ALP=ALPHA<  MT YP ) 

IF ( NSUCT <  EQ  <  1 .AND . DXX • GT . DGWT ) ALP=1 <  0 
IF<NSUCT.EQ.l)BPRES=BPREl*ALP*<l.>2.*AK0<MTYP>>/3'. 

IF ( NBPRES . EQ . 3 ) GO  TO  20 

IF(DXX.LT.0.01)G0  TO  30 

BL=BLEN 

BW=BWID 

BPR=BPRES 

IF<MRECT<EQ<DG0  TO  40  - 

BL=BLEN/2 . 

BU=BUID/2 . 

40  VE2=(BL**2r+BW**2.+DXX**2. >/<DXX**2. >  “ 

VE=VE2**0.5 
AN=BL*BU/(DXX**2.  ) 

AN2«=AN**2. 

ENM« <  2 . *AN*VE/ < VE2+AN2 ) )*< VE2+1 , >/VE2 
FNM*2.*AN*YE/<VE2-AN2> 

IF  <  MRECT . EQ . 1 ) BPR=BPRES/4 . 

AB=ATAN(FNM) 

IF<FNM<LT<0< )AB»PII+AB 
P < I > «P < I ) +BPR* < ENM+AB ) /P I I 


A8 


4750 

60  TO  70 

4700 

20 

DB*DXX/BUID 

4770“ 

* 

PS— .157-.22BDB 

4780 

IF (HRECT.EQ.O. AND.DB.LT. 2. 5 )PS=-.28*DB 

4790 

PS=10.**PS 

4800 

.  .  ... 

P(I)*P<I) 4BPRES4PS 

4810 

60  TO  70 

4820 

30 

P< I >  *P< I )+BPRES 

TOO' 

70 

DXX=DXX+DX 

4840 

10 

CONTINUE 

4890 

RETURN 

E 
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